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In the Standard Model (SM) neutrinos are left-handed 
massless particles and lepton number is conserved. How- 
ever, the strong evidence for neutrino oscillations [l[ indi- 
cates that neutrinos do have non-zero masses. An impor- 
tant question then arises regarding the origin of neutrino 
masses: whether they are of Dirac or Majorana type. If 
neutrinos are purely of Dirac type, they must have right- 
handed singlet components in addition to the left-handed 
states required in order to accommodate neutrino masses. 
In this case, lepton number is conserved. On the other 
hand, if there are Majorana-type neutrino states, a neu- 
trino cannot be distinguished from its own antiparticle. 
As a result, lepton-number-violating processes can occur 
in which lepton number changes by two units (AL = 2). 

There have been many experimental attempts to 
search for AL = 2 processes. The most thoroughly tested 
of these processes are neutrinoless nuclear double beta 
decays {01/(3(3) [2[. While the experiments are very sen- 
sitive, uncertainties in the nuclear matrix elements for 
0iy(3(3 would make it difficult to extract the mass scale 
of the neutrinos involved in such decays. As an alterna- 
tive, several authors have considered AL = 2 processes 
in meson decays [3|-l5| . 

The only existing experimental result for AL ~ 2 B 



meson decays is that of the CLEO collaboration, which 
searched for B+ — >■ h~£^£'~^ [6[, where h stands for tt, K, 
p, or K* and i stands for e or fi. They set upper limits 
on branching fractions for these decays in the range of 
(1.0-8.3) X 10-6 at 90% confidence level (CL) 0. Since 
6 — > c decays are in general favored in comparison to 
charmless B decays, it is interesting to extend the search 
to B~^ — >■ £)-£+£'+ decays. Two well-known diagrams for 
such decays are shown in Fig. [l](a) and (b). According to 
theoretical calculations, with a heavy Majorana neutrino 
of mass within the (2 — 4) GeV/c^ range, the branching 
fractions of _B+ -^ D^l^l''^ can be larger than 10"^ [j, 
|5[ with the diagram in Fig. [T] (b) giving the dominant 
contribution. 

In this paper, we report the first searches for the 
B+ -^ D-e+e+, D-e+fi+ and D' y+ tx+ decays. The 
results are based on a data sample containing 772 x 10^ 
BB pairs collected at the T(4S') resonance with the Belle 
detector at the KEKB [8|] asymmetric-energy e+e" col- 
lider (3.5 on 8 GeV). The Belle detector is a large-solid- 
angle magnetic spectrometer consisting of a silicon ver- 
tex detector, a 50-layer central drift chamber (CDC), an 
array of aerogel threshold Cherenkov counters (ACC), 
a timc-of-flight scintillation counter (TOF), and an ar- 



ray of CsI(Tl) crystals for an electromagnetic calorimeter 
(ECL) located inside a superconducting solenoid coil that 
provides a 1.5 T magnetic field. An iron flux-return lo- 
cated outside the solenoid is equipped with resistive plate 
chambers to identify muons as well as K^ mesons (KLM) . 
The Belle detector is described in detail elsewhere 19|. 

The analysis procedure is established using Monte 
Carlo (MC) simulations [10| , as well as data control sam- 
ples wherever possible. Since we have no prior knowl- 
edge nor widely-accepted model for the decay dynamics 
of i?+ — > I?~i?+£'+, the signal MC samples are generated 
uniformly over the three-body phase space, and we re- 
strict our analysis and interpretation to this model only. 

To reconstruct i?+ — >■ D~i^£'^ decays, we first look 
for an energetic same-sign dilepton and combine it with 
a D candidate requiring a proper charge combination for 
the dilepton. All charged tracks are required to originate 
near the interaction point and have impact parameters 
within 5 cm along the beam direction and within 1 cm 
in the transverse plane to the beam direction. 

Electrons are identified using the energy and shower 
profile in the ECL, the light yield in the ACC (A^p.e.) and 
the specific ionization energy loss in the CDC (dE/dx). 
This information is used to form an electron (£e) and 
non-electron (£e) likelihood. The likelihoods are utilized 
in the form of a likelihood ratio TZe = £e/{Ce + Ce) 



11| . Applying a requirement on TZe, we select elec- 
trons with an efficiency and a misidentification rate of 
approximately 90% and 0.1%, respectively, in the kine- 
matic region of interest. Muons are distinguished from 
other charged tracks by their ranges and their hit profiles 
in the KLM. This information is utilized in a likelihood 



ratio approach [1^ similar to the one used for the elec- 
tron identification (ID). We select muons with an effi- 
ciency and a misidentification rate of approximately 90% 
and 1%, respectively, in the kinematic region of inter- 
est. The efficiencies for electron (muon) ID are evaluated 
from data using the e~'"e~(/i "*"/!") pair production via the 



two-photon reaction 77 



(/i+/i ). Since the lep- 



ton ID performance is worse for lower-momentum tracks, 
we require the lepton momentum in the laboratory frame 
to be greater than 0.5 GeV/c and 0.8 GeV/c for electrons 
and muons, respectively. 

We require a same-sign lepton pair that has a total 
energy in the T(45) center-of-mass (CM) frame greater 
than 1.3 GeV. More than 95% of events have only one 
same-sign lepton pair. When there is more than one 
same-sign lepton pair, we choose the most energetic 
same-sign lepton pair from the three most energetic lep- 
tons in the event. 

Candidate D^ mesons are reconstructed in the D^ -^ 
K~^Tr~TT~ decay. Kaons and pious are selected from 
charged particles by applying hadron ID [13[. The 
hadron ID utilizes the time of flight measured in the 
TOP as well as iVp.c. and dE/dx in a likelihood ratio 
approach, which is similar to that used for lepton ID. 



We discriminate kaons (pions) from pions (kaons) with 
an efficiency of approximately 91% (95%) and a misiden- 
tification rate below 4% (6%) in the kinematic region of 
interest. The rates are evaluated from data using kine- 
matically reconstructed D*~^ — > D^tt"*" — >■ K~TT~^Tr'^ de- 
cays. The three tracks from the D^ candidate are fit to a 
common vertex and are required to have a K~^'!t~tt~ in- 
variant mass {Mktttt) within approximately ±10 MeV/c^ 
from the nominal D~ mass 14|. The Mk-ktt distribution 
is fit to two Gaussian functions with a common mean. 
The Mk-ht; mass window is chosen to be ±3 times the 
width of the narrower Gaussian component. The aver- 
age multiplicity of D~ candidates is 1.3 per event. If 
there are multiple D~ candidates, we choose the one with 
Mktttt closest to the nominal D mass. 

The same-sign dilepton and the D^ candidates are 
combined to form a B candidate, and are fit to a com- 
mon vertex. The B candidates are kinematically iden- 
tified using two variables: the energy difference, IS.E = 
Eb — i^bcam, and the beam-energy-constrained B meson 

is the beam en- 



mass, Mbc = V-^boam-Ps- Here, Sboam 

ergy and Eb and pb are the energy and momentum, re- 
spectively, oia, B candidate; these variables are defined in 
the CM frame. We select events with Mbc > 5.2 GcV/c^ 
and lAi^l < 0.3 GeV ("analysis region"). The signal re- 
gion is defined as 5.27 GeV/c^ < Afbc < 5.29 GcV/c^ and 
-0.055 (-0.035) GeV < A£; < 0.035 GeV for the e+e+ 
and e+jU+ modes (/x+/x+ mode), respectively. For back- 
ground studies, we use a subset of the analysis region 
that excludes the signal region ("background region"). 

One of the major backgrounds comes from the contin- 
uum production of quark pairs e+e^ ~> 99 (? = u, d, s 
and c) . The continuum background is discriminated from 
the signal by utilizing the difference of the event shapes 
in the CM frame. Since B mesons arc produced from the 
T(4S') resonance nearly at rest in the CM frame their 
final state particles are distributed isotropically. In the 
continuum, on the other hand, qq pairs hadronize back- 
to-back and give rise to a two-jet-like shape. To quantify 
the event shape characteristics, we use Fox- Wolfram mo- 
ments I15II with modifications optimized for exclusive B 
decays [l6|. A single discrimination variable, T. is ob- 
tained by applying a linear Fisher discriminant [17| to the 
moments and maximizing their discrimination power. 

In addition to J^, we also use the cosine of the polar 
angle of the B candidate fiight direction evaluated in the 
CM frame {cos 6b)- Since the T(4S') is a vector particle 
that decays to a pair of spinless B mesons, the cos 9b dis- 
tribution of the B mesons follows a jYnp oc 1 — cos^ 9b 
distribution, while random track combinations in the 
continuum have a nearly uniform distribution. 

The other major background comes from semileptonic 
B decays such as B ^ D't^v^X with D" — > X+tt^tt^, 
where X denotes any particle. Such decays can be mis- 
reconstructed as signal by combining a same-sign lep- 
ton from the decay products of the other B. In such 



background events, each lepton is produced along with 
a neutrino, resulting in large missing energy, while the 
signal tends to have small missing energy because there 
are no neutrinos in the final state. Here the missing 
energy, -Bmiss, is defined as Emiss = 2£;boain - Z^-^dot, 
where J2 ^dot denotes the sum of energies of all the de- 
tected particles in the event. Moreover, the same-sign 
leptons in such background events originate from differ- 
ent B mesons. As a result, the difference between the 
impact parameters of the two leptons in the beam direc- 
tion, Sz, tends to be larger in such background events 
than in the signal. Therefore, we use -Emiss and 5z as 
variables to suppress these backgrounds. 

The four variables, J-, cos 0b, E^iss and 5z, are 
combined together into a single likelihood ratio TZs ~ 
Ls/{Cs + £b), where £s(b) denotes the signal (back- 
ground) likelihood defined as the product of the signal 
(background) probability densities for each of the four 
variables. The two major backgrounds can be suppressed 
by applying a requirement on 7?,s. The probability den- 
sity functions (PDFs) are taken from the distributions in 
the MC samples. The background sample includes con- 
tinuum and BB components, where B decays are limited 
to & — > c decays. The optimal requirement on 7?.s is deter- 
mined by maximizing the figure of merit, eg/V-^b, where 
Eg is the signal efficiency estimated with the signal MC 
sample, and A^b is the number of expected background 
events in the signal region. Since only a small number 
of events remain in the signal region after the 7?,s re- 
quirement, the value of N\^ is obtained by scaling the 
number of events in the analysis region using the back- 
ground MC sample, where the scale factor is determined 
from the same MC sample but without the 7?.s require- 
ment. The optimal requirements on TZs eliminate more 
than 99% of the background while retaining 11-26% of 
the signal depending on the mode. 

In addition to the two dominant backgrounds de- 
scribed above, we checked backgrounds that might pro- 
duce a signal-like enhancement in the M^c-^E dis- 
tribution having more than one particle misidentified. 
Possible peaking backgrounds include i?+ — )- J/?/'(— > 
^+^~)i5r+7r+7r^, with the i~ and 7r+ misidentified as a tt^ 
and ^+ , respectively. Contributions from these decays are 
investigated using the MC sample that is approximately 
equivalent to 50 times the luminosity of the data sam- 
ple. The contribution of B^ — >■ D~h^h''^ decays with 
both same-sign hadrons {h^'') misidentified as leptons is 
estimated from the number of i?+ -^ D^h^h'^ events 
weighted by the h' misidentification rates, both evalu- 
ated in data. Background events from misreconstructed 
D~ mesons are studied using the D~ mass sideband. 
We studied charmless hadronic B meson decays as well 
as semileptonic B — >• Xu£iy decays using dedicated high- 
statistics MC samples, which are approximately equiva- 
lent to 21 and 14 times the luminosity of the data sample, 
respectively. 



After applying the TZs requirements, 5, 23 and 40 
events remain in the background region for the e"'"e+, 
e~^ fi~^ and fJ-^ fJ.^ modes, respectively. The background 
levels arc in good agreement with the expectations from 
the background MC samples; 4, 22 and 38 events, re- 
spectively. The signal region of the data sample is not 
examined until all the selection criteria are fixed and the 
systematic uncertainties are evaluated. From the MC 
samples the signal efficiencies are evaluated to be 1.2% 
- 1.9%, depending on the mode. Here the small differ- 
ence between the MC and data samples on the particle 
ID performance is corrected. In each case, the correc- 
tion is approximately 2% or smaller. The expected num- 
bers of background events in the signal region (N^^^) are 
0.18, 0.83 and 1.44 events for the e+e+, e+/^+ and fi'^n'^ 
modes, respectively. These background expectations are 
obtained by scaling the results of a two-dimensional fit 
to the background region, where we use a common back- 
ground shape for the three signal modes to compensate 
for the low statistics. The PDFs to fit the background 
distribution are an ARGUS function 181 for Mbr and a 



linear function for AE. We take the ratio of the integral 
of the PDF in the signal region to that in the background 
region; its value and error are 0.036 and 0.006, respec- 
tively. 

Figure [5] shows the Mbc- A_E distributions of events in 
the analysis region of the data sample, which pass all 
the selection criteria. The signal region is unblinded and 
no events are observed in any mode, which is consistent 
with the background expectations. Table H] summarizes 
the signal efficiency, the number of observed events and 
the expected number of background events in the signal 
region for each mode. 

The systematic uncertainties on N^^^ are also listed in 
Table m Each of the uncertainties combines the errors on 
the number of events in the background region and on the 
scale factor. For the latter each PDF shape parameter 
is varied by its fit error, and the resulting changes of the 
scale factor are added in quadrature. The fit procedure 
and the uncertainty evaluation are also applied to the 
background MC sample. Moreover, a mode- dependent 
PDF shape, taken from the background MC sample of 
each mode, is examined in the same manner. As a con- 
servative evaluation, the uncertainties obtained with two 
MC-bascd PDFs are added in quadrature in the uncer- 
tainty for each mode listed in Table ID 

Systematic uncertainties for efficiency determination 
are summarized in Table [Hi They are dominated by the 
tracking efficiency and the requirement of TZs- The un- 
certainty on the tracking efficiency is obtained by com- 
paring partially and fully reconstructed D*+ — > 7r+D°, 
D" — > Kg{-^ 7r+7r~)7r+7r~ decays in data and MC sim- 
ulation. The systematic uncertainties on the particle ID 
efficiencies are evaluated using the data control samples 
mentioned earlier. The uncertainty on the selection effi- 
ciency of the TZs requirements is evaluated from the ra- 



tio of the number of events in the signal region before 
and after applying the TZg requirement for data and MC 
samples using the B^ — > J/ijjK*^ mode. The number 
of events in the control sample is extracted by apply- 
ing the 2-dimensional fit described earlier with a PDF 
component for the corresponding decay. Since this con- 
trol sample does not represent the e~^ fi~^ mode very well, 
we take the larger of the two dilepton mode uncertain- 
ties for the e+/i+ mode. The same control sample is 
used to evaluate the uncertainty on the efficiency of the 
signal region acceptance. The same evaluation is ap- 
plied for the uncertainty on the efficiency of the Mktttv 
acceptance. A difference between the Mx-k-k shapes in 
data and MC would result in the different event frac- 
tions in the signal region. The control sample used is 
B^ — >■ £)"(—>■ K^'K~'K~)'n'^ , which is kinematically re- 
constructed after applying hadron ID requirements. 

No events are observed in the signal region. We set 
upper limits on the branching fractions based on a fre- 
quentist approach [l9j. We calculate the 90% C.L. up- 
per limit on the branching fractions including system- 
atic uncertainty, using the POLE program without con- 
ditioning [20[. Except for the uncertainty on iV^xp; ^-^ 
the systematic uncertainties, including those on the num- 
ber of BB events (Nn r ) and on the branching fraction 
of D^ — > i^+7r~7r~ [IJI, are assigned to multiplicative 
quantities in the upper limit calculation. These are found 



to be 8. 



and 9.7% for the e+e+, e^/i"*" and /i^/i^ 



modes, respectively, as summarized in Table HIl The 90% 
CL upper limits are (1.0 — 2.6) x 10^® depending on the 
mode, as listed in Table HI 

In summary, we have searched the Icpton-numbcr- 
violating _B+ -^ D^£+£'+ decays for the first time. We 
find no signal candidates. Assuming uniform three-body 
phase space distributions, we set the following upper 
limits on the branching fractions at 90% CL: B{B^ ^■ 
D-e+e+) < 2.6x 10"^, B{B+ -^ D-e+ii+) < 1.8x10"^, 
and B{B+ -^ D-fi+fi+) < 1.0 x 10^6^ 
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TABLE I: Results of the B+ -^ D~£+£'+ search; e is the signal 
reconstruction efficiency, Aobs is the number of events in the 
signal region, N^^p is the expected number of background 
events in the signal region, and U.L. is the 90% CL upper 
limit on the branching fraction. The efficiencies shown in the 
table do not include the branching fraction of the D^ decay. 



Mode 


e[%] 


A'obs 


Afbkg 
^ ^cxp 


U.L. [10-"] 


B+ ^ D-e+e+ 


1.2 





0.18±0.13 


< 2.6 


B+ -> D'e+^i+ 


1.3 





0.83±0.29 


< 1.8 


B+ -> D-fi+fi+ 


1.9 





1.44±0.43 


< 1.0 



[1] Y. Fukuda et al. (Super-Kamiokande Collaboration), 
Phys. Rev. Lett. 81, 1562 (1998); Y. Fukuda et al. 



(a) 



s+ 



(b) 




FIG. 1: Feynman diagrams for B+ -> D' £+£'+. 



TABLE II: Summary of multiplicative systematic uncertain- 
ties. The units are in percent. 



Source 


D-e+e+ 


D-e+fi+ 


D-fi+tj.+ 


MC statistics 


< 0.1 


< 0.1 


< 0.1 


Tracking efficiency 


5.2 


5.2 


5.2 


Lepton ID 


3.1 


3.5 


3.6 


Hadron ID 


1.4 


1.4 


1.4 


Mbc and AE 


2.0 


2.0 


1.5 


Mktttt 


2.4 


2.5 


2.4 


7^s 


3.0 


4.9 


4.9 


Nbb 


1.4 


1.4 


1.4 


B(D- ^ K+TT-TV-) 


4.3 


4.3 


4.3 


Sum 


8.8 


9.8 


9.7 
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FIG. 2: The Afbc-AE' distributions of D"e+e+(top), 
D~e^ fj,'^ (middle) and D^ fi^ jj.'^ (bottom) final states in data. 
The boxes indicate the signal regions. 



